Galaxy as dissipative N — body system. 
Chemical aspect. 
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Abstract. The evolving galaxy is considered as a system of baryonic fragments em- 
bedded into the static dark nonbaryonic (DH) and baryonic (BH) halo and subjected 
to gravitational and viscous interactions. Though the chemical evolution of each frag- 
ment is treated in the frame of one - zone model with instantaneous recycling, the 
local gas density and, therefore, its star formation (SF) activity is influenced by other 
interacting fragments. In spite of its simplicity the model provides a realistic descrip- 
tion of the process of galaxy formation and evolution over the Hubble timescale. 



1 The model 

The Plummer - type (PLT) Q halo is composed of DH and BH components. Their 
masses and softening parameters are correspondently equal to 10 12 Mq, 10 11 M Q , 
and 25 kpc, 15 kpc. The clumped baryonic matter (future galaxy disk and bulge) of 
mass 10 11 Mq is equally distributed among N = 2109 PLT particles with softening 
parameters 1 kpc. Initially particles are smoothly placed inside the sphere of radius 
50 kpc, which rigidly rotate around z axies with initial angular velocity that provide 
rotational support of matter in its outer equatorial region. They are involved into the 
Hubble flow (Hq = 75 km/ s-Mpc~ l ) and have random velocities 0-^10 km/ sec. Their 
subsequent dynamics is defined by particle - halo and particle - particle gravitational 
interactions and by viscous BH influence and is followed by means of effective N 

- body integrator with individual time step. Gravitational terms are defined in a 
standard manner. The comparison of model results with data of more sophisticated 
SPH approaches shows that viscous acceleration term can be choosed as AY ISC = 

- Vi • | V, | /R VISC ■ P BH(r)/pvisc ■ mr/imr + mf ar ). Here p BH is the 
baryonic halo local density and Vi is a particle velocity vector, pvisc = 0.2 cm -3 
and Rvisc — 10 kpc are model parameters. A particle mass is defined as a sum 
of gas mf as and stellar mf ar components which depend on its evolutionary status. 
Initially mf tar = 0. To account for the influence of baryonic halo rotation onto 
the dynamics of baryonic particles in the disk galaxy plane A VISC is set equal to 
(0.15,0.15,1.0)- A VISC . 

The characteristic time step in the integration procedure for a particle is Sti — 
10~ 2 • min^-y/j R,; 3 | / | A i3 |, | Ry | / | V y - |], where the V i3 = Vi — Vj and the 
Ay = Ai — Aj is a difference of accelerations of i and j particles. Such step value 
provides a nice momentum and energy conservation over the integration interval of 
about 15 • 10 9 year. 

The local gas density pi is defined in a SPH - like manner as a ratio of mass of 
gas contained in the sphere of radius Hi around the i th particle to its volume. The 



Hi is defined (using the quicksort algorithm) by the requirement that the number of 
particles inside this sphere is equal to Nb — 21 (« 1 % of N). The particle eligible to 
star formation (SF) events is choosed as particle which still have a sufficient amount of 
gas (exceeding 0.01 of particle initial mass) and densities exceeding the critical value 
PminSF = 0.01 cm -3 during the fixed time interval At$F = 50 Myear. Such active 
particle produces stars with a Scalo - type initial mass function with SF efficiency e is 
defined as e = 0.5 -pi/10 cmT z ■ (1 — ex<p{— 10 cm -3 / pif). Each SF event is followed by 
SN explosions which inject into the particle the enriched gas of metallicity Zi + AZ. 
Initially Zi = 0.0 and AZ = 0.01 as the average value of enrichment for stars having 
metallicities in the range Z = 0.001 -f- 0.04. The returned fraction of gas is set equal 
to R = 0.25. This gas is assumed to be immediately mixed with old (heavy element 
deficient) gas. 



2 Conclusion 

The rapidly rotating protogalaxy formes a three - component system resembling a 
typical spiral galaxy: a thin disk and spheroidal component made of gas and stars 
and a dark matter halo. The disk component poseses a typical spiral galaxy rotation 
curve and the distribution of radial and z - th velocities of particles clearly shows the 
presence of the central bulge. The metallicities and the global metallicity gradient 
resemble distributions observed in spiral galaxies. During the first 2 Gyr of evolution 
only about 20% of protogalaxy baryonic component is transformed into the stars, 
but, finally, the fraction of the stellar population grows up to 92%. The galaxy global 
star formation rate (SFR) can be considered as a continuous succesion of short bursts 
which initially doesn't exceed 28 M Q yr^ 1 . Then SFR gradually decreases up to the 
value of about 2 Mq yr~ x . All these data as well as surface densities distributions of 
stellar and gaseous components are in nice agreement with present date observational 
data of our own Galaxy ||, || and with results of more sophisticated theoretical 
approaches [@, ||, [| |l| . 
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